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Contributions to the bound-state dynamics of fermions in local quantum field theory from the 
region of large relative momenta of the constituent particles, are studied and compared in two 
different approaches. The first approach is conventionally developed in terms of bare fermions, 
a Tamm-Dancoff truncation on the particle number, and a momentum-space cutoff that requires 
counterterms in the Fock-space Hamiltonian. The second approach to the same theory deals with 
bound states of effective fermions, the latter being derived from a suitable renormalization group 
procedure. An example of two-fermion bound states in Yukawa theory, quantized in the light-front 
form of dynamics, is discussed in detail. The large-momentum region leads to a buildup of over- 
lapping divergences in the bare Tamm-Dancoff approach, while the effective two-fermion dynamics 
is little influenced by the large-momentum region. This is illustrated by numerical estimates of the 
large-momentum contributions for coupling constants on the order of between 0.01 and 1, which is 
relevant for quarks. 
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I. INTRODUCTION 



The notion of a bound state of fermions is based mainly 
on the examples of atoms and nuclei. The common fea- 
ture of these systems is that they are non-relativistic. 
This means three things. Kinetic energies of the fermions 
are small in comparison to their rest mass energy. Dom- 
inant interactions are not able to produce fast moving 
fermions from the slow ones and hence no significant 
large-momentum spin effects are generated. Creation 
of additional particles can be neglected and one can de- 
scribe the bound states as built from a fixed number of 
fermion constituents. These features are all related to 
the fact that the domain of large relative momenta is not 
important, which is easiest to describe in mathematical 
terms in the case of bound states of two fermions, such as 
positronium or deuteron. Their wave functions are self- 
consistent solutions to the non-relativistic Schrodinger 
equation H\9) = ^l^") with Hamiltonian H = Hq + Hj, 
where Hq denotes the kinetic energy operator and Hi 
stands for the interaction operator. The matrix element 
(12|_ff/|l'2') is the quantum Coulomb, or Yukawa po- 
tential with a repulsive core, respectively. The ket |12) 
denotes a state of two fermions labeled 1 and 2, with 
all their quantum numbers collectively denoted by these 
labels. The self-consistency of this well known picture 
means that the wave function ^'(1,2) = (12|5') quickly 
vanishes when the relative momentum of fermions in the 
bound-state rest frame of reference, p = pi — p2, ap- 
proaches values comparable with the fermions' reduced 
mass /i (the speed of light, c — 1). 
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The success of the Schrodinger picture for bound states 
of fermions extends also to quarks. This is reflected in 
the constituent quark model (CQM) being the primary 
means for classification of hadrons in the particle data 
tables, and providing a benchmark for more advanced 
approaches. However, the self-consistency of the non- 
relativistic picture is much harder to maintain for bound 
states of quarks up, down, and strange, than for systems 
such as positronium and deuteron. This is because the 
hadronic wave functions tend to have considerable com- 
ponents with \p\ comparable to or exceeding fi for the 
quarks, and the domain of large relative momenta be- 
gins to play a significant role in the binding mechanism. 
One is also interested in the description of hadrons mov- 
ing with speeds very close to the speed of light. Since 
the fast moving hadrons and their interactions cannot 
be consistently described within a purely non-relativistic 
Schrodinger framework, theorists use the Feynman par- 
ton model in that case. Unfortunately, the binding mech- 
anism of partons is far beyond the scope of the parton 
model. As alternative to these models, one can approach 
the issue of bound states of fermions using quantum field 
theory (QFT), where the corresponding operator H ap- 
pears to contain all the relevant information about rel- 
ativistic effects in the domain of large relative momenta 
of the constituents. 

However, the relativistic description of bound states of 
fermions in QFT, especially in quantum chromodynamics 
(QCD) in the case of light quarks, makes the conceptual 
difficulties with the constituent picture even greater than 
in the simple models. In the equation Hqcd\'^) — ^-l^); 
all factors remain unknown. Major reasons for this sta- 
tus of the theory originate in the large relative momen- 
tum domain in the motion of virtual particles. This is 
illustrated by the following examples. One is that lo- 
cal QFTs lead to canonical interaction Hamiltonians Hj 
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that change individual bare particle energies by unlimited 
amounts (spin-dependent factors grow with momentum 
transfers). The large momentum range is enhanced and 
leads to divergences, invalidating the concept of a non- 
relativistic picture entirely unless special conditions, such 
as an extraordinarily small coupling, are met. Another 
example is that the interactions create new bare parti- 
cles and this effect contributes to the boosting of bound 
states, which implies that the motion of bound states 
is associated with multi-particle components and the ad 
hoc limitation to a fixed number of bare constituents goes 
out the window in relativistic QFTs. A third example 
is that even the state with no constituent particles, i.e. 
the ground state of a theory, or vacuum, proves to be 
so complicated that no approximate solution of verifiable 
accuracy has been conceived yet, although many ansatzes 
have to be and are employed in practical attempts. In 
these circumstances, the main theoretical approach to 
bound states of quarks (mainly heavy ones that move 
slowly) is based on the lattice version of QCD, and great 
progress has been achieved in numerical studies of the 
theory that way. Nevertheless, it appears that a quan- 
titative explanation of how the constituent picture with 
a simple Hamiltonian could be an approximate solution, 
remains a conceptual and quantitative mystery. No de- 
tailed constituent wave function picture for relativistic 
field quanta in Minkovsky space has been theoretically 
identified or is expected to readily follow from the lattice 
approach alone. The question of convergence of the bind- 
ing mechanisms in the domain of large relative momenta 
of constituent particles, remains open. 

This article provides some numerical arguments that 
the required constituent picture with well-controlled 
large relative momentum domain may become in prin- 
ciple identifiable if one provides a precise definition of 
the constituents as effective particles, in distinction from 
the bare quanta of the local theory. Thus, the process 
of solving a theory is arraiiged in two steps, which is 
typical in lattice approach [n], or sum rules 0. In the 
first step, one derives an effective dynamics, and, in the 
second step, one attempts to solve the effective theory 
instead of dealing directly with the original degrees of 
freedom. Here, one derives the effective fermions of size 
using a suitable boost-invariant perturbative renor- 
malization group procedure for their Hamiltonians. The 
procedure is carried out up to second-order perturbation 
theory, and the resulting dynamics is compared with the 
standard picture, where the finite scale A is absent. In 
distinction from the diverging bare dynamics, the effec- 
tive one comes out limited to the momentum range given 
by A, and this scale is reduced using differential equa- 
tions to the most suitable one for description of bound 
state properties in terms of a fixed number of the corre- 
sponding constituents. In the renormalized Hamiltonian 
picture, the point-like bare particles of the local theory 
correspond to A = oo and their dynamics heavily involves 
large momenta, and multi-particle states, for any finite 
value of the coupling constant. However, the situation is 



completely changed when A is lowered to values compara- 
ble to the bound state masses. The binding is described 
by a new Schrodinger equation, Hx\'^) = ^^l^*), where 
the Hamiltonian, Hx — Hqx + Hjx, is written in terms of 
creation and annihilation operators for the effective par- 
ticles, 6^ and bx for fermions, and dx for anti-fermions 
(and a\ and a a for bosons). This effective particle picture 
is discussed in the present paper. 

The key physical reasons for the hope that the effec- 
tive constituent picture does emerge from QFT can be 
understood by recalling again what happens in the well- 
known cases of atoms (or positronium) and nuclei (or 
deuteron) discussed above. These systems can be under- 
stood in terms of constituents for quite different reasons. 
The explanation of the difference is limited below to the 
positronium and deuteron, but the two examples are suf- 
ficient to make the point that concerns all bound sys- 
tems of fermions in QFT. In the Schrodinger quantum- 
mechanical picture of positronium, the coupling constant 
that appears in the Coulomb potential is very small in 
comparison to 1, i.e. e^/An = a ~ 1/137. Therefore, the 
interaction produces quite small binding energy, of or- 
der and the e'^e~ bound-state mass is dominated 
by 2me. The relative- motion wave function is propor- 
tional to (a^/i^ -I- |p p)~^, independently of the fermion 
spins. When one extends this picture by embedding it 
in quantum electrodynamics (QED), one finds out that 
the initial wave function is so small for large momenta p, 
that no significant correction is able to emerge from that 
region and alter the original picture with the Coulomb 
potential. This is found by expanding the theory term 
by term in powers of a around the initial Schrodinger pic- 
ture. The interaction linear in a (Coulomb force) is suf- 
ficient to describe the main features of fermionic bound 
states in QED, and higher powers of a are not important 
for theoretical understanding of the bulk of the bound 
state structure. Although the integrals in the corrections 
run over the momentum range that formally extends far 
beyond /z, the coupling constant is too small for the rel- 
ativistic fermion spin factors and particle creation pro- 
cesses to produce any major modifications of the leading 
approximation. This feature will be farther discussed in 
the next Sections. 

In the meson-exchange models of the deuteron binding 
mechanism, the analogous coupling constant is three or- 
ders of magnitude larger than in QED. If one attempted 
to use QFT to derive the Yukawa potential using the 
same strategy as in QED, and to calculate corrections, 
the perturbative procedure would fail. The interactions 
would accelerate nucleons to the speed of light almost im- 
mediately on the bound-state formation time-scale, many 
new particles would be created, and the large momentum 
relativistic "corrections" would dominate the "leading" 
non-relativistic terms. One could then ask why the fa- 
mous one-boson-exchange (OBE) potentials, such as the 
Yukawa potential with a core, could still be used in phe- 
nomenology of relativistic nuclear physics and work self- 
consistently in the nuclear bound-state equations any- 
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way. What saves the picture of a fixed number of rel- 
atively slow nucleons interacting through exchange of 
mesons from serious inconsistency when one includes the 
elements of QFT, is that the interactions responsible for 
emission and absorption of mesons by nucleons contain 
form factors that limit momentum transfers to values so 
small that the interactions are effectively weak, much 
weaker than a change of a in QED by the factor 1000 
would imply. Consequently, the binding energy is much 
smaller than the sum of two nucleon masses, e.g. about 
2.2 MeV for deuteron. The wave functions of such gen- 
eralized relativistic nuclear physics picture are not over- 
whelmingly extending into the large-relative-momentum 
domain because the form factors eliminate coupling to 
that region and the non-relativistic Yukawa potential 
with a repulsive core is not invalidated by huge correc- 
tions. It could not be so with bare point-like fermions in 
local QFT, but it does work in the phenomenological pic- 
ture of effective particles with the form factors. By the 
way, this example is not intended to suggest that nucleon 
dynamics should be completely derivable directly from a 
local QFT that ignores the existence of quarks. In fact, 
a scenario for how to derive the effective nuclear physics 
picture from QCD will be discussed in the last Section 
of this work. Nevertheless, the nuclear physics picture 
does indicate that an effective particle dynamics may in- 
volve large coupling constants in potentials that resemble 
perturbative second-order interactions with form factors. 

In QCD, none of the two schemes can apply separately. 
On the one hand, the effective coupling constant in the 
constituent QCD picture cannot be as small as in QED, 
because QCD is characterized by asymptotic freedom, or 
infrared slavery. This means that the effective coupling 
strength is expected to grow when the scale of relevant 
momentum transfers decreases. The coupling constant is 
already on the order of 0.1 at transfers on the order of 
100 GeV and it may be much larger for transfers on the 
order of nucleon masses. Therefore, the effects that have 
marginal size in the eigenvalue problem for approximate 
Hqed, such as spin splittings, are expected to be much 
larger and more important for understanding eigenstates 
of Hqcd, and the initial approximation is not as sim- 
ple as in QED. On the other hand, one cannot freely 
insert form factors into the local Lagrangian for quark 
and gluon fields because it would spoil the local gauge 
symmetry structure. The contact with QCD would be 
irreversibly lost. 

The situation is changed when QFT is approached us- 
ing the idea of similarity renormalization group proce- 
dure for Hamiltonians especially in the case of QCD 
Q , and when one combines the similarity idea with the 
concept of form factors in the Hamiltonian interaction 
vertices for effective particles ||]. Initially, the coupling 
constant is small due to asymptotic freedom and one can 
think of using the small coupling constant in canonical 
QFT to solve the renormalization group equations for 
H\ using a perturbative expansion. The method avoids 
small energy denominators in the perturbative calcula- 



tion entirely and the non-perturbative part of the dynam- 
ics remains untouched in the perturbative calculation of 
H\. That way one derives effective- particle Hamiltoni- 
ans that involve vertex form factors of small width A in 
the interaction terms. One can have sizable couplings in 
H\ with small A, as required by infrared slavery, without 
losing control over the size of corrections to the leading 
constituent picture in diagonalizing H\ . The spectrum of 
such Hamiltonians can be sought numerically because the 
form factors limit the range of momenta strongly enough 
for a discrete computer code to cover the pertinent re- 
gion, as in the nuclear physics case discussed above. This 
idea has already been studied qualitatively in a simple 
numerical matrix model |^ using Wegner's flow equation 
(and its generalizations). The more detailed effective 
particle calculus used in the present work in the case of 
Yukawa theory, is already known to produce asymptotic 
freedom in H\ for QCD Q . The new approach has been 
also extended to the whole Poincare algebra ||^ in QFT. 

The present paper provides elementary numerical es- 
timates of the orders of magnitude of the interactions 
that appear in QFT in the bound state dynamics of two 
effective fermions of size A^^. Yukawa theory is used to 
avoid complications related to gauge symmetry while one 
still preserves some of the singular large-momentum com- 
ponents in the spinor factors that characterize fermions. 
The well-known issue of triviality in Yukawa theory is 
irrelevant here since our goal is to estimate the size of 
corrections in the bound-state dynamics in an effective 
theory, rather than in the ultra-violet (UV) dynamics of 
the initial QFT. The Yukawa example serves only as a 
source of typical UV factors that QFTs provide anyway, 
no matter if the theory is trivial, asymptotically free, or 
otherwise. 

The key qualitative question is by how much the 
Hamiltonian H\ derived in QFT might differ from 
familiar models, especially from the non-relativistic 
Schrodinger picture with a Yukawa potential (or a 
Coulomb potential in the case of exchange of massless 
mesons), for given values of a and A. Another question 
is related to the fact that the exact solution of renor- 
malization group equations for H\ and subsequent exact 
diagonalization of H\ should lead to spectra that are in- 
dependent of A. However, when one calculates H\ in 
perturbation theory of low order, such as the second or- 
der that characterizes the Coulomb and Yukawa poten- 
tials, the A-dependence must appear. Bound-state en- 
ergies may depend on A when A is made too small or 
a is made too large. The question is how large is the 
residual A dependence of second-order corrections to the 
Coulomb-like picture in how large a range of values of a 
and A that can be self-consistently (i.e., without signif- 
icant A-dependence) reached in lowest orders of pertur- 
bation theory, for Hx itself in the renormalization group 
part of the calculation, and for the eigenvalues and wave 
functions in the bound state perturbation theory around 
the Coulomb- like ansatz. Both questions are answered 
in the next Sections by estimating the size of those cor- 
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rections that are most important for large momenta, and 
which would lead to divergences in the absence of A. The 
results imply that the most dangerous corrections that 
might diverge in the absence of A, turn out to be quite 
small even for sizable coupling constants. This provides 
some initial ground for the hope that the effective par- 
ticle dynamics can be attempted numerically in QCD in 
the same way. An outline of how the same approach, if 
successful, could then in principle also lead from effec- 
tive HxQCD to the effective meson-exchange interactions 
between nucleons in nuclear physics, is also provided. 

Many farther dynamical issues of great interest, con- 
cerning effective fermions and their binding in QFT, can 
be posed and investigated using the effective particle 
calculus, but none are studied here in detail. For ex- 
ample, the case of phonon exchange between electrons 
as approximated by a QFT is not discussed here [0, 
and symmetry restoration issues in not-explicitly sym- 
metric Hamiltonian approach through coupling coher- 
ence is also not considered. Instead, this work is 
focused on the readily accessible quantitative estimates 
that show the magnitude of the difference between the 
convergent two-effective fermions bound-state dynamics 
and the diverging dynamics of two bare fermions in the 
approaches based on the Tamm-Dancoff truncation in lo- 
cal QFT llj, In fact, this paper starts with the 
Tamm-Dancoff- like approach to local theory, because it 
is more familiar. Then, the effective particle approach 
is introduced, with its new features exposed through 
the contrast with the Tamm-Dancoff one. Thus, Sec- 
tion II provides definitions of the rcnormalized Tamm- 
Dancoff scheme in Subsection II A, and the effective par- 
ticle scheme in Subsection [I B . Both approaches in- 



volve a universal procedure for obtaining a two-fermion 
eigenvalue equation that can be compared with the non- 
relativistic Schrodinger equation. This universal proce- 
dure is for brevity called reduction. I t is denoted by the 
symbol R and described in subsection [I C. Section |l| in- 
troduces a bare light-front Hamiltonian in Yukawa quan- 
tum field theory that serves as a starting point for sub- 
sequent Sections. The canonical Hamiltonian is supplied 
with some regularization factors ta and counterterms. 
Section IV describes details of the approach that treats 
bound states of fermions as if they could be viewed as 
made of two bare fermions. This approach is shortly 
called approach 1 and encounters large-momentum con- 
ceptual and calculational difficulties that are removed 
by switching over to the approach discussed in detail 
in Section Namely, Section ^ presents the effective 
fermion approach, called approach 2, using comparison 
and contrast with the approach I. In the approach 2, 
bound states of fermions are treated as built from ef- 
fective fermions of size A^^. Conclusions are drawn in 



Section VI. Appendixes contain information about the 
operation R, spinor factors, momentum variables, and a 
Coulomb bound state wave function used in the calcula- 
tions. 



II. DEFINITIONS 

This Section provides definitions of two light-front 
Hamiltonian approaches to the bound state dynamics, 
the rcnormalized Tamm-Dancoff approach (approach 1), 
and effective particle approach (approach 2). They are 
compared through numerical estimates in next Sections 
in the case of a bound state of two fermions. The pro- 
cedure of reducing a theory to the two-fermion sector 
is mathematically the same in both cases, except for a 
different definition of what is meant by the two-fermion 
states. In the approach 1, the starting point will be iden- 
tical with bare fermions in QFT, in the approach 2 it will 
be the effective ones. The reduction procedure, denoted 
by _R, is described at the end of this Section for both 
approaches simultaneously. 



A. Renormalized Tamm-DancofF approach 



This approach is represented by the following diagram. 



C 



Hr 



Hi 



solve 



(1) 



0) The initial step on the left denotes a canonical 
derivation of a field theory Hamiltonian from its La- 
grangian, quantization, and normal ordering with respect 
to the bare vacuum state, dropping all diverging terms 
on the basis of hindsight knowledge that the normal or- 
dered Hamiltonian will eventually contain counterterms 
of the same structure, but well defined. 

1) The next step is called regularization, which is neces- 
sary because the canonical Hamiltonian leads to ultravio- 
let divergences due to the infinite range of energy scales to 
sum over in the intermediate states when one calculates 
physical observables (such as scattering cross-section) us- 
ing perturbative formulae. One imposes an ultraviolet 
cutoff (we denote it by A) on the range of momenta that 
are included in the sums. Observables calculated with 
such regularized Hamiltonian explicitly depend on both 
the parameter A, and on the way the regularization is im- 
posed. For example, the regularization may preserve or 
violate some symmetries - an oblate cutoff function would 
not be rotationally symmetric, or a frame-dependent cut- 
off would not respect Lorentz symmetry, and finite cor- 
rections due to such regularizations would be required. 
To remove the artificial dependence of observables on reg- 
ularization, one has to add new terms to the Hamiltonian 
(called counterterms and denoted X^) that also depend 
on the regularization. The construction of counterterms 
is based on the idea of recovering the contribution from 
dynamics above the cutoff, which is to lead to finite and 
cutoff independent results when A ^ oo. The coun- 
terterms structure is thus determined by the structure 
of divergences. The latter can be computed and X^ is 
than adjusted to cancel the cutoff dependence. It is less 
straightforward to remove finite effects of regularization, 
but this issue will not be important in the next Sections. 
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The regularized Hamiltonian : -ff^,„ : +X is denoted 
by H^. 

ii) The last arrow indicates solving of the eigenvalue 
equation for H^. A two-steps procedure is used. 

Step a) First one finds eigenstates of whose domi- 
nant component for vanishingly small coupling constants 
is equal to one bare fermion. These states represent what 
one could call a physical fermion. The solution is found 
from the eigenvalue equation for the whole by re- 
ducing that equation with the help of operation R to 
an equivalent equation for the Fock component with one 
bare fermion. If one requires the eigenvalue to be finite, 
one has to include in a mass counterterm of a calcu- 
lable form. 

Step b) Then, one makes a reduction R of the eigen- 
value problem for to a two-bare-fermion subspace, to 
find an eigenstate of the Hamiltonian that is domi- 
nated by a pair of bare fermions for infinitesimally small 
coupling constants. The parameters in the resulting two- 
fcrmion eigenvalue problem are expressed in terms of the 
physical fermion mass found in step a) above. It turns 
out that the calculated eigenvalues still depend on the 
cutoff (some diverge if A cx)), although the individual 
matrix elements of the reduced two-body Hamiltonian do 
not depend on the cutoffs once one includes mass coun- 
terterms calculated in ii a). Therefore, there is a problem 
of how to construct countcrterms that would remove A 
dependence from physical results [p^ . 

It is important to mention here an alternative ap- 
proach, which is also inspired by the Tamm-Dancoff trun- 
cation but includes renormalization group idea in a more 
sensible way Namely, one can assume that some 

effective fermion representation of the eigenstates of 
does exist. One tries then to build a triangle of renormal- 
ization (the notion of the triangle in the case of similarity 
renormalization group scheme for light-front Hamiltoni- 
ans can be found in Ref. for a sequence of calculations 
with growing cutoffs A and growing numbers of fermions 
and bosons, adjusting countertcrms to the momentum 
cutoffs differently in different Fock sectors. The hope is 
to understand regularities in the obtained structures that 
produce stable outputs when the momentum cutoffs and 
the particle numbers are made so large in the triangle 
that one can follow the renormalization group flow down 
from the large values and recognize interactions that have 
a universal character. Moving then along the renormal- 
ized Hamiltonian trajectory in the triangle, one would 
be able to arrive at the desired effective particle picture. 
One of the key tools one would use in studying the tri- 
angle would be the coupling coherence [Tl| | that requires 
that the interactions retain their forms when one moves 
around the triangle with large values of the cutoffs. This 
approach is thus relying on regularities that need to be 
discovered in a vast amount of intricate analytic and nu- 
merical data. The data can also be interpreted with the 
help of perturbative techniques wherever possible. As a 
variation of this idea, one may even attempt to include 
new degrees of freedom, such as Pauli-Villars pseudo- 



particles, and treat them as a suitable way of writing 
down required countertcrms |p^ . However, the latter 
strategy faces a problem of sending the regulator masses 
to infinity (it is not clear how to do it using computers) 
or living with the additional degrees of freedom that are 
absent in the initial theory. One also uses discretization 
of the longitudinal light-front momenta (usually denoted 
by P^) |17| , which limits the number of constituents. 
However, the transverse dynamics is not smoothed out 
that way and poses the same divergence problems. The 
initially helpful longitudinal cutoff has to be relaxed in 
order to recover boost invariance. At that stage, the 
same problems reappear fully again and remain unsolved. 
Thus, the already known options for general searches for 
effective fermion dynamics are hard to pursue. In con- 
trast, in the present article, it is taken for granted that 
one can assume from the very beginning a large degree 
of rigidity in the derivartion of effective theory by writ- 
ing it in terms of creation and annihilation operators for 
effective particles. The assumption allows one to narrow 
the search for effective dynamics using a unitary renor- 
malization group procedure , which is the basis for the 
approach 2 introduced below. 



B. Renormalized effective particle approach 

The procedure consists of three steps that are de- 
scribed later in detail in Section and can be repre- 
sented by the following diagram. 

£ : ff,,„ : ^ : H^n : ^ Ha ^ solve (2) 

The steps 0) and i) are the same as before except that 
one works with the bare creation and annihilation oper- 
ators for efficient book-keeping for Hamiltonian terms at 
all times, instead of storing a huge number of selected 
matrix elements of H^. Note that many other operators 
could have the same selected matrix elements. Also, ma- 
trix elements of operators change when the basis states 
change even if the operators remain unchanged, and the 
same operator may have different matrix elements when 
different states are used to represent constituent parti- 
cles. 

ii) This step, marked with U\ in the diagram, is made 
using the procedure of renormalization group for effec- 
tive particles (RGEP) ||, ||. The procedure is defined 
as a unitary rotation of creation and annihilation opera- 
tors by an operator U\. Hamiltonians H\ are expressed 
in terms of the effective-particle creation and annihila- 
tion operators that depend on the "width" parameter A. 
A ranges from cxi in to a finite value on the order of 
bound-state masses in the effective constituent dynam- 
ics. That such small scales can in principle be reached 
using perturbation theory |6| , is possible because RGEP 
procedure is designed according to the rules of the simi- 
larity renormalization group for Hamiltonians . It sep- 
arates the perturbative from non-perturbative aspects of 
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the theory (see the original articles). By the same design, 
the Hamiltonian H\ cannot change invariant masses of 
effective-particle Fock states by more than about A in a 
single interaction. Thus, the emission of effective bosons 



FIG. 1: The main difference between and H\ is that 
the latter is for effective particles of size A~^, instead of the 
point-hke bare ones that the former is for, and interaction 
terms in H\ are limited by form factors fx. The form factors 
are marked in figures by a black blob. 



by effective fermions is possible only if the associated ki- 
netic energy change of relative motion of the particles 
does not exceed A. Consequently, when A is small, Fock 
sectors with different numbers of effective particles are 
coupled weakly even for sizable coupling constants, like 
in the nuclear physics example discussed in the Introduc- 
tion. 

in) This step is analogous to the step ii) in the ap- 
proach 1 and amounts to solving the eigenvalue problem 
for effective Hamiltonian H\. The key difference, how- 
ever, is that when one works using the basis of effec- 
tive particles in the Fock space, states with two effective 
fermions couple only to states with similar relative mo- 
menta. Therefore, the large relative momentum domain 
remains suppressed, and it can be described using pertur- 
bation theory without assuming that the coupling con- 
stant is very small. Thus, when one solves the eigenvalue 
equation for H\, one can introduce a new perturbation 
theory for the reduction operator i?, expanding in powers 
of Hi\. This gives an equivalent Hamiltonian that acts 
only in the dominant Fock space sectors. There are two 
steps to do, as in the approach 1. 

In the step a), one first considers eigenstates dominated 
by one effective fermion, which defines a physical mass 
of a physical fermion in the approach 2. Then, in the 
step b), one finds an equation describing bound states 
of two effective particles. The parameter A is a key to 
the procedure. Its value determines whether derivation 
of the effective Hamiltonian H\ and its reduction by the 
operation i? to a model subspace Hamiltonian, denoted 
by Hfj, is possible in perturbation theory. The smaller 
is A the simpler are the approximate solutions for bound 
states of effective fermions, in the sense that they tend 
to reduce to the dominant effective Fock sector. But if 
A is too small, the step ii) of derivation of Hx in pertur- 
bation theory loses accuracy (the perturbative integra- 
tion of renormalization group equations begins to signif- 
icantly cut into the bound-state dynamics). Therefore, 
A cannot be lowered too far using perturbation theory 
for H),. The optimal choice of A is the one that com- 
bines the simplest perturbative expansion for H\ with 



least complicated computer diagonalization of H\. The 
main criterion for choosing the right range for A is that 
the calculated observables are not sensitive to variation 
of A over that range. 

The final comment concerns Refs. 18, where a 
different concept of a bound state calculus has been de- 
veloped using coupling coherence in second-order pertur- 
bation theory for Hamiltonian matrix elements, also in 
the similarity scheme but without the constraint to a 
boost-invariant unitary rotation of creation and annihi- 
lation operators (from the bare to effective ones) . In dis- 
tinction from these works, the approach 2 is not based 
on the coupling coherence because no coherent structure 
is known a priori in the region of small A, far from the 
initial canonical structure. Instead, one uses a perturba- 
tive expansion for the effective-particle renormalization 
group flow in terms of a suitably defined coupling con- 
stant and tries to find out the relevant structures in a 
prescribed basis, in which the expansion in powers of 
the coupling constant may be extrapolable to its physi- 
cal values. Here, the analysis of Yukawa theory is limited 
to numerical estimates of the size of only a few impor- 
tant terms of first and second order in a bound-state per- 
turbation theory, which goes beyond the renormalization 
group flow. In the flow itself, it is already known that 
higher-order calculations are certainly possible even in 
much more complex theories (see e.g. ||]). The effective 
bound-state dynamics is virtually unknown ||20|. 



C. Reduction procedure 

The following scenario occurs several times in next Sec- 
tions. There is an eigenvalue equation for a Hamiltonian 
H = Hf) + Hi, 



(3) 



which is too large to solve exactly on a computer in the 
sense that the number of a priori important basis states 
is inflnite. One looks then for an equivalent Hamiltonian 
that acts only in a limited subspace of states. One way of 
constructing the model subspace dynamics is to use the 
transformation R ^^[^ Details of the transformation 
are written in Appendix lAI. The idea is following. 

One denotes the projection operator on the chosen sub- 
space of the whole Fock space F by P, and the projector 
on the complementary space, 1 — P, by Q. If the inter- 
action Hamiltonian Hi is small in the sense that it only 
weakly couples states from the subspace PF to states in 
the subspace QF, then one can calculate an operator R 
that produces certain eigenstates of H from eigenstates 
of a new Hamiltonian Hfi that has eigenstates contained 
in the subspace PF. The transformation R leads to the 
following expression for the Hamiltonian Hn acting in 
the subspace PF, expanded in powers of the interaction 
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Hamiltonian Hj (cf. Eq. 



III. CANONICAL THEORY 



= (^| [PHP + \PHI 



Q 



2 Ej — Hq 



Note that the Hamiltonian Hn does not depend on the 
eigenvalues of H, but only on the eigenvalues of Hq. 
Namely, -ffoN) = Ei\i). In particular, one can de- 
fine Hq in conjunction with the subspace PF so that 
Hi — H — Hq is as weak as one can get simultaneously 
with preserving control over the spectrum of Hq. For ex- 
ample, in the case of two fermions, it is often useful to 
take Hq equal to a sum of the kinetic energy and Coulomb 
potential operators when the fermions interact through 
emission and absorption of massless bosons with a small 
couphng strength. The strength of the coupling (size of 
Hi) may be small when the coupling constant itself is 
small, or when a form factor allows only small range of 
momentum transfers from the fermions to bosons and 
between the fermions. 



We consider a theory of fermions of two kinds inter- 
acting through a Yukawa coupling with scalar bosons. 
The common starting point for all approaches to a two- 
fermion bound-state dynamics discussed in this paper is 
the Lagrangian: 



(5) 



where ^5 = (-0'^-'^^ i/j^^^) . That is, C includes the kinetic 
terms for two kinds of fermions (both of mass m), the 
kinetic term for massless scalar bosons and a point-like 
(Yukawa) interaction of the fermions with bosons. 

The canonical Hamiltonian is defined using light-front 
quantization p3| , where the evolution of states in "time" 
~ x'^ + x'^ is generated by a three-dimensional integral 
of the energy- momentum-density-tensor component 
over the light-front hyper-plane a;"*" = 0. The Hamilto- 
nian is given by 



dx (fx^ 



(0 



r 



',(0 



(6) 



ipm are fermion fields that satisfy free Euler-Lagrange 
equations with mass m. Quantum commutation relations 
for all fields are satisfied by replacing all fields at a;+ = 0, 
[a;-*- = {x^ , x'^) and x~ = x'^ — x^ are the "spatial" coor- 
dinates] , by their Fourier expansions in terms of creation 
and annihilation operators (see Appendixes ^ and ^ for 
details). The canonical Hamiltonian derived this way 
leads to infinities in perturbative calculations of observ- 
ables, due to summation over intermediate states with an 
infinite range of relative momenta that are reached by the 
local interactions from any state with only finite relative 
momenta. In order to regulate the Hamiltonian en bloc, 
including the binding mechanism, one has to introduce 
a cutoff A on the relative momenta in the interaction 
terms. The cutoff requires counterterms that remove the 
cutoff-dependence from observables in the limit A ^ oo. 
Thus, the Yukawa-theory Hamiltonian takes the form 

H^ = Hq + H^ + H^ + X^ , (7) 

where Hq denotes the contribution of the first two terms 
in Eq. (0), Hy of the third term, and of the fourth 



one. X represents the counterterms. To be specific, 
the relative momenta (see Appendix ^ for definitions of 
X and K^) of created or annihilated particles are limited 
by inserting factors ta = exp {—n^'^/A^)rs[x) in the in- 
teractions terms, e.g. see Eq. (|^). In the current study 
of Yukawa theory, the x-regulator factor rs {x) is needed 
only for the mass-counterterm construction. In all other 
formulae, it can be replaced by 1. Since rs{x) does not 
appear in the final results, its form is left unspecified here 
(cf. §). 

The free part, Hq, is 

+ / b]^«^« ■ (8) 

There are three interaction terms. The one relevant here 
consists of emission and absorption of a boson by the 
fermions, i.e. 
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'^q"p2<y2 Pl^l P2<^2 ' 



9 P20'2 Plt^l '^Pl'^l 



/i.e. 



(9) 



The other two are: those that create fermion-antifermion 
pairs from a boson and vice versa, and the instantaneous 
interaction between fermions and bosons mediated by 
fermions. These two interactions are not written exphc- 
itly here since they do not contribute to equations ana- 
lyzed in the following sections. It is quite possible, and 
actually desired, that these additional interactions will 
improve the constituent picture when included in the dy- 
namics. However, it cannot happen in calculations of the 
lower order than 4th, which is not known yet. 



IV. APPROACH 1: BOUND STATES OF TWO 
BARE FERMIONS 

This Section reviews the renormalized Tamm-Dancoff 
procedure for two-fermion bound states. One starts with 
a single fermion eigenvalue problem, and then proceeds 
to the two-fermion case. 



A. One-fermion eigenstates 

The one-fermion eigenvalue equation is first obtained 
by assuming that the coupling constant in the theory is 
infinitesimally small and the dominant part of the eigen- 
state is provided by a single bare-fermion Fock state. The 
quantum numbers of the lowest-mass eigenstate corre- 
spond then by definition to one physical fermion asso- 
ciated with the fermion field in the initial Lagrangian. 
When the coupling constant is made finite and grows, the 
eigenvalue equation is not soluble exactly with currently 
known mathematical methods, and one has to investigate 
results that follow from various attempts to find approx- 
imate solutions. One such attempt is made by reducing 
a cut-off dynamics to the one-bare fermion Fock sector 
for finite coupling constants, too. In that c ase, the pro- 
jection operator in the operator R (see Sec. p C] and Ap- 
pendix has the form P = E<t /W4!^^|0>(0|4!^ For 
finite cutoffs and sufficiently small coupling constants g, 
one uses expansion in powers of g to evaluate the corre- 
sponding Hn. Up to order g^, this leads to an equation 
Hnlk) = P-\k), with 



,-L2 



P- 



'A 



k+ 



(10) 



where 



results from emission and re-absorption of 



bosons and Xjf is contributed by the counterterm pro- 
portional to wh. Since rrA is a diverging function of A, 



Xff has to be adjusted to remove that effect. Note that 
should not and does not depend on the fermion mo- 
mentum components k'^ and fc-*- in the light-front the- 
ory developed here. The expression for the eigenvalue 
P^ is unique to light-front Hamiltonian dynamics with 
regulators preserving boost symmetries, so that for all 
momenta k-^ and k'^ there is one and the same value of 
vrif. Common approaches to field theories in standard 
time quantization lead to momentum-dependent "mass 
terms", i.e. the candidates for mf in the eigenvalue 
Ef{k) = (mj + |fcp)^/^ depend on k and boost symme- 
try is violated. It is also common in those approaches to 
consider only slowly moving fermions, especially at rest 
with respect to the frame distinguished by the definition 
of H'^, and to identify the value of the function mf{k) 
at with the desired fermion mass. Such approach is 
at best only approximate in the case of fc ^ and cre- 
ates conceptual difficulties in description of bound-state 
constituents with large relative momenta. Therefore, the 
standard quantization is replaced here by the light-front 
scheme. 

The result ( [lO| ) for requires a counterterm of the 
form 



A2 
2 



— / dx-r'i 



A2 

4m^ log ^ -I- C 



,(11) 



where the constant C is a finite part of dimension m^. 
This condition removes A-dependence from the physical 
fermion mass my in the limit A — > oo. 

B. Two-fermion bound states 



Reduction (Sec p^I C] ) to a two-bare fermion Fock sector 
employs (note that the two fermions are selected to be of 
different kinds) 

p=Y: j[PW2]h^M^'^'\ombM\ (12) 

and leads to an equation for an eigenstate of Hn^ which 
can be written as 



\Pcj) = / biP2]P+(2^)^<53(P-pi -p2)</'.6Wt5(2)t|o). 

(13) 
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The eigenvalue equation can be then written in terms 
of the two-body wave function (p^ = 4>a{k) (cr denotes 
spin quantum numbers of both fermions a = {ai,a2\, 



where Ek = {k^ + m^)^/^, and the potential kernel vqbe 
will be discussed below. 

Note the mass m/ in Eq. (p^, i.e. the physical fermion 
mass obtained from the earlier reduction to one bare 
fermion space, Eq. (|lO|). Expanding the bare mass m, 
in both the integration measure factor {EiE2)~^^'^ and 
potential vqbe, in a series of powers of g around nif, 
leads to an equation in which only physical mass to/ 
enters. This way the bound-state dynamics for two bare 
fermions is related to a physical fermion mass parameter. 
This step makes a connection between the bare fermions 
in the two-body problem and a physical fermion obtained 
in the one-body reduction in the previous subsection. 

To simplify farther discussion of the two-fermion eigen- 
value equation, we denote the single-fermion mass- 
eigenvalue TO/ by TO (i.e. the subscript / is dropped). 
The two-fermion bound-state mass Mfuii can be rewrit- 
ten as Mfuu = 2m — Eb- When Eb <C to, the eigenvalue 
takes the form (M|„„ - 4m2)/4m = -Eb + E%/Am « 
—Eb- Therefore, the eigenvalue on the right-hand side 
of Eq. ( p^ ) can be thought of as the binding energy Eb ■ 

Since the regulator function ta respects kinematical 
boost invariance of the light-front scheme, this equation 
is independent of the total momentum of two fermions. 
There is also no explicit A-dependence in the matrix el- 
ements of the potential vqbe in the limit A — s- cx): 

liniA^oo 1'OS_E(o-lO-2(T3(T4fcfc') = 

2to2 q+ \ba bcjp^g^^ 
+the same|Fig.|f, , (15) 

where a = q^/Att. The notation adopted here is taken 
from Ref. [|| (see also caption of Fig. || here). The spinor 
matrix elements are, 

UiU2 = , x\[m{xi + X2) - a^cr^ixik^ - X2ki)]x2, 

^XiX2 

(16) 

where = [Ij 0] or [0, 1], depending on the fermion-spin 
projection on the z-axis. 

The potential (|l^) is a quite complicated, nonlocal 
function of fermion momenta. But in the region with 
both momenta fc, k' <C to, it simplifies to the well-known 



for definition of fca that forms k together with fc-*-, see 
Appendix ^ . Namely, 



(14) 
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(a) (b) 

FIG. 2: Two kinds of terms in one boson exchange poten- 
tials. Following Ref.Q, the initial (rightmost), intermediate, 
and final (left-most), states are denoted by c, b, and a, re- 
spectively. For example, in diagram (a), pjj^ = , p^^ = 
p+, ba = M|+5 - m\ be = M|+5 - m^ q+ = p+ . 



Coulomb potential (cf. Appendix P), 

1 ^ ^ 

VCoulomb = -4:TTa-^ — ^— . (17) 

[k — k'Y 

Unfortunately, this heuristic result is not meaningful be- 
cause the region of large relative momenta of the fermions 
introduces important corrections. From Eq. (|l^) one can 
see that for one of the relative momenta (fc or k') much 
bigger than the other, and than the fermion mass, the 
spinor factors become proportional to the larger of the 
two momenta. For example, if k' ^ fc,TO, one obtains 
uu ^ k' , and two such factors compensate the denom- 
inator that grows as k'^ [Q. The potential becomes a 
function of xi and X2, being a constant in the transverse 
momentum directions. A constant potential in the trans- 
verse momentum space is a two-dimensional 5-function 
potential in configuration space. Such potentials with a 
negative coefficient lead to bound states of infinite nega- 
tive energies in the non-relativistic Schrodinger equation, 
and the light-front transverse dynamics is of this type. 
One could try to rely on the regulators ta with finite 
A to remove the problem, which would correspond to 
smearing of the 5-potential in position space. The eigen- 
values of the equation would then depend on A. One 
could naturally try to make the coupling constant g de- 
pend on A. However, the interaction is specific to the 
Fock sector under consideration, it is much more compli- 
cated than a (5-function itself, at least by the presence of 
the additional a;-dependent factors, and it is unlikely that 
a change of 5 to a function of A can remove the cutoff 
dependence from all eigenvalues. The procedure cannot 



-^^W+E(43 \ ^^V0BBi^y,k,k')^Ak')= ^ Mk) 
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be based on exact solutions in seeking simple remedies 
for the cutoff-dependence. 

A perturbative search for suitable counterterms in the 
two-body sector was carried out up to second order in 
Ref. iQ. The problem of buildup of overlapping di- 
vergences in the two-body sector was in principle solved 
for low-mass eigenvalues in the limit A ^ oo in Ref. 
p6| . However, that procedure was not able to deal with 
the small-energy denominators in a perturbative calcu- 
lus for effective Hamiltonians. It is worth stressing that 
the cutoff-dependence problem appears not only when 
one sends A to infinity, but already for finite A, since 
the eigenvalues have interpretation of and some of 
them become negative and lose physical interpretation 
when A and the coupling constant are large enough. This 
problem concerns lowest angular momentum states. Note 
also, that in some of the states, divergences may cancel 
out p5[ | due to interference of various spin amplitudes in 
the limit A oo. Nevertheless, the two-body equation 
as it stands is not convergent in the large relative mo- 
mentum domain, and the cutoff-dependence invalidates 
the non-relativistic approximation as a means for seek- 
ing a conceptually satisfying solution of the divergence 
problem, especially for sizable coupling constants. 

A calculation described in the remaining part of this 
Section illustrates how the overlapping divergence prob- 
lem arises in approach 1 in a quantitative way. An anal- 
ogous calculation with effective fermions in Subsection 
VB2, will not produce such diverging results. This dis- 
tinction is important because the effective-particle ap- 
proach allows then a self-consistent farther development 
of the Fock-space expansion and inclusion of higher-order 
terms in the Hamiltonians, which together may eventu- 
ally lead to a well-founded constituent picture. 

Since the potential in the (non-relativistic) region of k 
and k' small in comparison to m has the Coulombic form, 
one can ask with what accuracy Eq. (|lj) can be ap- 
proximated by a Schrodinger equation with the Coulomb 
potential (given in Appendix y). The potential can be 
re-written in the form 



VOBE = VCoulomb + 



(18) 



and corrections induced by Av estimated in perturba- 
tion theory. The Coulomb potential does not depend 
on spins of the interacting fermions. Therefore, in 0-th 
order of the bound-state perturbation theory, there are 
four degenerate states with the lowest mass M, and the 
same momentum-space wave functions: a triplet of spin-1 
states, and a singlet of spin 0. 

To estimate the first-order energy correction one has 
to find eigenvalues of the 4x4 matrix of matrix elements 
{ipoilAVltpaj) 7 where i and j refer to the different spin 
configurations. The eigenstates of this matrix have spin 
structure: {U + it), (Ti - it), TT and ij. The lowest 
mass eigenstate is (ti — it)- The fist-order correction to 
the Coulomb energy for this state varies between numbers 
of the order of -4 IQ-^Eo for a = 0.01 to -0.09£:o for 
a = 0.6. Note that a is also present in the wave function 



(j)Q and the results are not connected by a straightforward 
multiplication by the ratio of the coupling constants, al- 
though both results are small, indeed. For a greater then 
0.6 the first-order correction would continue to be rela- 
tively small, but the second-order corrections (see below) 
become unacceptably large, and this is why we do not 
discuss a significantly larger then 0.6. 

Unfortunately, in the second order a convergence prob- 
lem in the domain of large relative momenta of fermions 
destroys self-consistency of the naive perturbative pro- 
cedure around a non-relativistic approximation. To see 
this in a transparent way, one can make a number of sim- 
plifications and isolate the origin of corrections that grow 
with A, without worrying about details of secondary im- 
portance. The point is that the Coulomb basis functions 
have quickly falling off tails in momentum space. A tail is 
still small but greatly enhanced by first-order corrections, 
and the second-order correction involves matrix elements 
that already diverge with the cutoff A. 

To see the origin of the overlapping large-relative- 
momentum divergence in the second-order energy correc- 
tion, one needs to analyze matrix elements of the type 



AS(2) = 



\Av- 



En — Hn 



Vco 



-Av\(j)o 



(19) 



Such elements involve integration over four relative mo- 
menta of fermions: the leftmost wave function argument 
denoted by fc; , the momentum of states between the left 
Av and l/{Eo - Hq - Vcouiomb), denoted by pi, the mo- 
mentum between the operator l/{Eo — Hq — Vcouiomb) 
and right Av, denoted by Pr, and the argument of the 
right wave function, denoted by kr- The matrix element 
can be split into a sum of 2^ parts, with each part distin- 
guished by saying how large is each of the four integrated 
momenta in comparison to the fermion mass m, smaller 
or larger. 

Firstly, since the Coulomb wave functions strongly 
limit their arguments, a part with ki and kr large is a 
very small contribution in comparison to the part with fc; 
and kr small. Therefore, one looks for important contri- 
butions assuming that ki and kr lie within several widths 
of the Coulomb wave functions. An ad hoc number used 
here is Aam. 

Secondly, there is no bound-state wave function lim- 
iting the intermediate momenta pi and Pr, and integrals 
over them extend up to the cutoff A ^ oo. Eq. ( |l^ ) 
shows that for large momenta the spin-flip part of the po- 
tential dominates other parts. This dominant part is se- 
lected here and denoted by Au||, both fermions have op- 
posite spin orientation and both have their spins flipped 
in the interaction. For the purpose of estimating the 
order of magnitude of the large-momentum spin-flip con- 
tribution, pi and Pr are considered larger than m and 
{EQ-Ha-Vcouiomby^ is replaced by - I /Ha, neglecting 
terms that would vanish when a — *■ 0. Then, the resol- 
vent becomes diagonal in momentum space and pi — Pr 
are commonly denoted by /c2 ■ Details of how the cutoff 
A was initially introduced are not important for the or- 
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der of magnitude estimate. Therefore, the cutoff function 
is shghtly changed to simplify the integration. Namely, 
the initial ta limits changes of invariant masses in each 
of the vertices, see Fig. y, producing a complex shape 
of the /c2-integration boundary with details that depend 
on the small momenta fc; and fc^, irrelevant to the di- 
vergence issue at hand. The main role of the cutoff in 
Ai5^^^ is to provide the upper limit on the range of inte- 
gration over k2. Therefore, one can estimate the size of 
the large-momentum range contribution by introducing 
a new cutoff kmax, equal to the maximum value that fc2 
can take, and let kmax — ^ oo when A ^ oo. The depen- 
dence on kmax will indicate dependence on A. By the 



way, introducing a kmax in Ai?^^) is also a way one could 
introduce cutoffs if regularization were first imposed in 
the reduced eigenvalue equation ( |T^ ) rather than in the 
initial QFT Hamiltonian, cf. Using the simplifica- 

tion of TA to 9 (kmax ~ I ^2 1 ) here should not be understood 
as advocation of that procedure in a precise QFT calcu- 
lation. 



As a result of these steps, the large relative-momentum 
part of the second-order energy correction can be esti- 
mated from the following expression (the tilde indicates 
the simplifications made in AE^^^^ as discussed above). 



J 



4am 



(Pki_ 

7K 



d^k2 Z"*"™ d^k 



^Mki)Avti{ki,k2)—Avu{k2,k3)Mh) (20) 

E3 -Ho 



r 



The range of integration over k2 in this expression, is 
shown in Fig. ^. Since the potential Aw|| approaches 



4a m 



k m 4a m 

MAX 

FIG. 3: The shaded area represents the range of integra- 
tion over I fell and Ifel in Eq.(|20[) with the potential matrix 
Aw||(fci,fc2). In the lower-right corner, both momenta are 
equal zero. 




FIG. 4: Dependence of the most singular part of the 2nd- 
order correction (Eq. po| ) on the cutoff kmax- This correction 
diverges logarithmically for large kmax even for small coupling 
constants a, though the matrix elements of the two-body 
Hamiltonian do not depend on kmax for momenta smaller 
than kmax • 



a constant for k2 ^ m, one can expect a logarithmic 
dependence of Ai?^^^ on kmax, 

d'^k2 ^ - 1 ^ ^ kmax 

-—Av-^l{ki,k2)—Av^l{k2,k3) - log . 
ii2 -Ho m 

(21) 

A numerical evaluation of the 12-dimensional integral 
produces an estimate of the actual size of the logarithmi- 
cally diverging correction. The results for different values 
of the coupling constant are given in Figure I The er- 
rorbars indicate the standard deviation of a Monte Carlo 
routine used in the computation. All other parts of the 



second order two-fermion bound-state mass correction 
(i.e. the parts with external momenta bigger than Aam, 
or internal momenta smaller than m, or parts without 
change of the fermion spins) cannot compensate this di- 
vergence. Note that the corrections can quickly reach the 
order of 10% for coupling constants of the size expected 
in quark physics when the cutoffs are made larger than 
100 quark masses. Even if one assumes that to ~ 1/3 
GeV, momentum transfers on the order of 30 GeV would 
be already far too large to fit under such cutoffs. The 
constituent wave-function picture would encounter seri- 
ous conflicts with first principles of the underlying theory 
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much earlier. If one attempted an analysis starting from 
much smaller values of the quark masses, closer to the 
Standard Model, the allowed cutoffs could become ex- 
tremely small and the procedure would be stuck with 
huge corrections from the large-relative momentum re- 
gion, out of control. 

An analysis analogous to the approach 1, though more 
advanced with respect to the actual bound state dynam- 
ics (thanks to solving the singular eigenvalue equations 
numerically) was carried out in Ref. ||lj]. Those au- 
thors used the renormalization group ideas and calcu- 
lated in perturbation theory two-body counterterms that 
could remove the cutoff dependence from the bound-state 
eigenvalues to large extent. Their procedure assumed 
that eventually the renormalization triangle with grow- 
ing cutoffs and particle number could be obtained in a 
much larger project with many Fock components, but 
they did not go beyond their first step since the coun- 
terterms in higher orders (when the potential would be 
more complicated) would be much more difficult to find 
and understand. Also, if the counterterms are guessed 
sector by sector, a clear connection to the initial field the- 
ory is not directly visible and great determination would 
be required to pursue studies of such complex eigenvalue 
problems and their dependence on the cutoffs. It is also 
worth stressing that for a small enough, there exist eigen- 
states of the low-angular momentum two-body problem 
with finite eigenvalues even for A ^ oo This 
does not change the fact that important eigenvalues in 
the bound-state equations for bare constituents depend 
on the cutoffs when the latter are much larger than the 
fermion masses. It is hard to guess a self-consistent path 
out of the problem without systematic understanding of 
many relativistic QFT effects that currently appear to lie 
far beyond a few-body approximation. 

Nevertheless, one could try to farther develop stan- 
dard perturbative methods of constructing counterterms 
using Wilson's renormalization group idea [^2[ in the ap- 
proach 1. But in the presence of spin-induced divergences 
in the bound states of fermions, there exists a problem 
of degeneracy of states that, most probably, cannot be 
resolved without use of the newer idea called the sim- 
ilarity renormalization group procedure This new 
method avoids the small energy denominator problem in 
degenerate perturbation theory. It can be used in evalu- 
ating effective Hamiltonians with small renormalization 
group parameters, and these parameters can then play 
the role of cutoffs that are sufficiently small to define the 
effective constituent dynamics self-consistently. Since the 
similarity idea is also one of the key ingredients of the 
effective-particle approach used in the next Section, it 
is important to point out here in general terms how it 
differs from the standard procedure of integrating out 
high-energy degrees of freedom for bare particles. 

In the presence of divergences, the standard require- 
ment is that counterterms are chosen in such a way that 
there is no A-dependence in a model Hamiltonian acting 
in a subspace of limited "kinetic" energy, see Fig. |^. The 



(a) 



(b) 




FIG. 5; (a) Tlie standard renormalization procedure is based 
on reduction of a space of states. This leads to small energy 
denominators in perturbation theory, (b) Similarity renor- 
malization procedure is based on a rotation of basis and avoids 
small denominators in perturbative derivation of effective dy- 
namics, because it integrates out only these energy changes 
that are larger than A. 



problem is that if there is no finite energy gap in the basis 
of eigenstates of Hq between the retained and integrated 
out states, then small energy differences between states in 
the retained subspace and outside of it appear in denom- 
inators in perturbative transformation R that is repeat- 
edly used in the renormalization group part of the whole 
calculation. The small denominators lead to large num- 
bers multiplying terms which depend on A, and intruder 
states can spoil the utility of that procedure completely. 
Therefore, to cancel the A-dependence by the calculated 
counterterms one would have to make incredibly precise 
calculations for very large A. In the similarity approach, 
the construction of counterterms and the transformation 
to small cutoff Hamiltonians are defined quite differently, 
see Fig. ||b. Instead of integrating states out, one brings 
the Hamiltonian matrix elements close to the diagonal 
in the sense that the new interactions change eigenval- 
ues of Hq only by less than an arbitrary but finite A. 
The small energy denominators are avoided because one 
never eliminates transitions with energy changes smaller 
than A, and so the denominators are always bigger than 
A. When one comes to the step of reduction R of H\ to 
a few-body bound-state eigenvalue problem, one already 
has a well defined effective theory with a small width A 
and, in principle, no trace of the diverging cutoff param- 
eter A. 

One more advantage of the similarity approach is that 
one can also use it to introduce the notion of effective 
particles and switch over to the eigenvalue problem 
written in terms of them, instead of the bare ones. This 
option allows us to perform practical calculations in QFT 
in a new way described in the next Section. 
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V. APPROACH 2: BOUND STATES OF TWO 
EFFECTIVE FERMIONS 

This Section briefly reviews the RGEP for deriving 
Hamiltonians Hx for effective particles of size and 
then applies H\ in the Yukawa theory to a bound state of 
two effective fermions. The presentation refers to some 
steps made in the previous Section, but in the context 
of effective particles instead of the bare ones. Key dif- 
ferences are pointed out, that lead to the two-effective 
fermion dynamics that converges in the region of large 
relative momenta. 



A. Renormalization Group for Effective Particles 

The RGEP is defined through a unitary rotation for 
creation and annihilation operators |9 



(22) 



Each operator (such as a Hamiltonian) can be expressed 
in terms of both sets of the operators: b or b\, and it has 
different matrix elements in the Fock-space basis built 
using the operators of each kind. The idea of the RGEP 
is to perform the rotation (p^ ) in such a way, that the 
Hamiltonian expressed in terms of b\, called the "effec- 
tive Hamiltonian of width A", denoted H\, contains ver- 
tex form factors fx of width A in all interaction terms. 
There are infinitely many interaction terms in Hx for all 
finite values of A and their strengths vary with A Q . The 
choice for fx made here is. 



/a = exp 



(M, 



created annihilated / 



A4 



(23) 



tive expansion for Hx , on top of the error margin result- 
ing from approximations used to solve the Schrodinger 
equation with Hx- On the one hand, one tries to get down 
to as small A as possible so that the non-perturbative di- 
agonalization will require smallest possible range of en- 
ergy scales to handle explicitly, using a computer. On the 
other hand, one expects errors due to use of perturbation 
theory in evaluating Hx to grow with reducing A, and A 
should not be too small. The reason is that if A ^ 0, the 
Hamiltonian becomes almost diagonal, which is equiva- 
lent to solving the non-perturbative dynamics of bound 
states, and a perturbative calculus for Hx must fail at 
some point before A becomes equal to the scale of the 
non-perturbative phenomena. 

The same rotation Ux provides also means for con- 
structing counterterms in H^ — Hx=oo- A Hamiltonian 
Hx with a finite A has a band diagonal matrix in the effec- 
tive basis, and each effective Fock basis state is directly 
coupled only to a limited set of other effective states that 
have energies within the range of A. Therefore, the ef- 
fective theory splits into a chain of theories that couple 
only to near neighbors, without jumping up to arbitrar- 
ily large scales such as A in the approach 1. Conse- 
quently, if counterterms introduced in H^ cause that no 
A-dependence appears in the matrix elements of Hx in 
the effective basis states, there can be no A-dependence 
generated in observables calculated in perturbation the- 
ory to any finite order when A/ A tends to infinity. This 
will be seen in detail at the end of this Section, but it can 
be observed already here that the bare-particle approach 
described in the previous Section did not have this prop- 
erty and A-dependence could show up in the bound-state 
dynamics even though matrix elements of the Hamilto- 
nian reduced to the two-body sector did not depend on 
A for all finite relative momenta. 



where Mcreated is a total free mass of all particles created 
by a given term in Hx, and M annihilated is a total free 
mass of particles annihilated by the term. 

If the unitary transformation Ux were known exactly, 
there would be no A dependence in the spectrum of 
Hx- But when Ux (and Hx) are calculated in pertur- 
bation theory, the approximation leads to some residual 
A-dependence of theoretical predictions for observables. 
The sensitivity of results to variation of A, provides a sim- 
plest test for how large errors one makes in the perturba- 



1. Effective Hamiltonian - 0th and 1st order 

The only change in the Oth-order Hamiltonian (free 
part, order g^) is that the bare operators such as b^'b are 
replaced by b\bx- In order g^, the effective Hamiltonian 
has the form, 



H 



A(l) 



2 . 

H H / [PiP2g] 2{2tt)^S^ {pi -p2-q) exp (-k^^^^/A^) rs{xg/xp, ) exp 



2—1 <7l ,(T2 



A4 



^ '"c(A^P2CT2A^piCTlA"P20-2'"pitTi 



+ h.c. 



(24) 



Note that expressing 6's by 6a 's induced the form factor factor that depends on A is equivalent to 1 when A/ A 
fx in Hx- This form factor causes that the regularization 
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oo. 



2. Effective Hamiltonian - 2nd order: mass term 



One can calculate the term in H\ of order 



9 



that 



contains bxh\ and see that it contains a mass-squared-like 
term with a divergent A-dependence. Therefore, one has 
to add a counterterm to the initial Hamiltonian that has 
exactly the same form (^Tj) as in the approach 1. After 
including this counterterm, the form of the effective mass 
term in H\ is (in the limit A — s- cx)) 



H\Sr, 



(25) 



Note that the renormalization is carried out now at the 
level of full theory in the whole Fock space, not after re- 
duction to a specific Fock sector. Therefore, for example, 
there are no sector-dependent mass counterterms. Since 
the regulators did not violate any kinematical light-front 
symmetries, the calculated mass term does not depend on 
particle momentum (i.e. the relativistic form of the dis- 
persion relation does not change, there is only a change 
of the value of the effective fermion mass). 



3. Effective Hamiltonian - 2nd order: potential term 



where 



5m\ 



167r2 



exp 



d.-(l 



2{z - m^) 



A4 



6m 



2\2-\ 



const. 



(26) 



Second-order terms in Hx that contain two creation 
and two annihilation operators for effective fermions do 
not contain any dependence on A when A — > cxi, and 
no counterterms are needed of such form. Therefore, the 
complete answer for these terms is 



(2) 

XWa (Pl'P2,P3,P4,Cri, (72, (73,(74) , (27) 

where 

^(pi,P2,P3,P4,o-i,o-2,(73,(74) = - g'^uiU2U3U4-^ fac^2{a, b, c)\Fig^ + the same|Fi(,.|& , (28) 

I 



and.F2(a,6,c) = [{xba+il~y)bc)/{ba^ + bc^)]{fbafbc-l), 
adopting conventions from Ref. |^ with x ~ Pi /{Pi + 
p^) and y = pt /{pt ~^Pt)- Despite that Fig. ||is referred 
to in order to define notation, just like in approach 1, the 
potential v\ is quite different from the OBE potential 
of Eq. (^5|). For example, note the different kind of 
denominators and the presence of the key formfactors 

fac- 



4- Effective Hamiltonian - 2nd order: other terms 

In the second-order effective interaction, there are also 
other terms, similar to the terms shown above, or de- 
scribing interactions with explicit participation of bosons. 
There are also terms creating or annihilating two addi- 
tional particles. None of these terms contribute to the 
second-order effective equation that will be obtained in 
the following Section in the case of bound states of two 
effective fermions. In general, the RGEP allows one to 
construct both the transformation U\ and Hamiltonian 



H\ in perturbation theory to an arbitrary order in g. Cal- 
culating higher-order corrections is ultimately the only 
way for finding out how large corrections they produce. 
The approach 2 is limited here to terms of the second 
order, because it was the case in the approach 1. 



B. Solving eigenvalue problem with H\ 

In the case of bound states of two effective fermions, 
the reduction procedure is based on the same rules as 
in the approach 1, except that the effective particles in- 
teract with vertex form factors of width A and the large 
relative-momentum convergence is improved. Also the 
change of particle number is severely limited in strength, 
since massive particles cease to be produced when A is 
lowered below their mass, and the emission of massless 
particles changes energies by amounts linear in the ex- 
changed momentum. The changes of order k are larger 
than the changes of order k x k/m in the fermions' en- 
ergies when A is smaller than m. The departure point of 
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the process of solving the bound state dynamics is the 
eigenvalue equation for the single fermion states. 



1. Reduction to one effective fermion subspace 



This step produces an equation Hjf \k) — P |fc), where 



k+ 



(29) 



and m?j is the physical fermion mass, by definition of the 
same value as in the approach 1. It comes out indepen- 
dent of A by the virtue of adjusting once and for all the 
mass squared counterterm in . The same adjustment 
involves fixing the free finite constant in Eq. (26) so 
that for some value of A = Aq the physical fermion mass 
eigenvalue to/ equals to the experimentally found num- 
ber. The interesting point is that the same eigenvalue is 
subsequently automatically obtained for all values of A 
and the physical dispersion relation satisfies all require- 
ments of special relativity. This is the simplest manifes- 
tation of the general rule that physical results should be 
independent of A, as it is only a parameter of a unitary 
rotation of the basis. 



2. Reduction to two-effective fermions 

Using transformation R to reduce Hx to the two- 
effective particle subspace without restrictions on the rel- 
ative momenta, one obtains a quantum mechanical inter- 
action that can change the invariant mass of the two par- 
ticles by a certain A if, and only if the interaction acts 
more than A/ A times. Thus, the approach 2 produces 
an effective Hamiltonian which is free from the overlap- 
ping divergence problem discussed in Ref . p6| , and in the 
previous Section in approach 1. However, in order to 
make a connection with the non-relativistic two-particle 
Schrodinger quantum mechanics, which was not system- 
atically available in the approach 1, one needs now to 
limit the relative momenta in the effective two-particle 
Fock sector to fc < z, where z is a new parameter re- 
quired for defining the new operation R that enables one 
to define the procedure of introducing the non-relativistic 
limit. 

Therefore, a new transformation R is now defined to 
lead to a model Hamiltonian Hn that acts only in the 
subspace of the two-effective particles Fock sector with 
limited invariant masses (Fig. |6|). Thus, not only the 
number of effective particles is limited, but also the range 
of their relative momenta. It is required that Hj^ has the 
same spectrum of low lying energy levels as Hx has in the 



whole space. This step is no longer related in any way 
with infinite renormalization problem as in the approach 
1. The existence of such reduction is plausible only be- 
cause Hx has a small width A and in subsequent orders 
of perturbation theory in Hi\, corrections to the effec- 
tive potential result from the coupling to an additional 
region of relative momenta, which is always limited by A 
in every new order. 




R 




EFFECTIVE QFT 



MODEL QM 



fa) 



fb) 



FIG. 6: Reduction of effective QFT, (a), to relativistic quan- 
tum mechanics in a model subspace, (b), is possible in the 
approach 2 thanks to the form factors in the effective inter- 
action vertices. 



The projection operator used here is 

P = E / biP2]6lf 6(f |O)(O|6(^>«0 [z ^ \k\) 



(30) 



where k is the relative momentum of effective particles of 
momenta pi and p2. Although introduction of z is use- 
ful from the conceptual point of view, the formfactors fx 
imply that z is not important in practice if only the low- 
est order (i.e. g^) model Hamiltonian is calculated (see 
Fig. H). I would, however, affect the model Hamiltonian 
in or der through terms such as the last two terms in 
Eq. ([A4| ). 

The effective Schrodinger equation has then the form of 
Eq. (Q) with vqb e replaced by a new potential, denoted 
vrxi which is a sum of two terms (Fig. ob). The first term 
is the projection of H^^j-^^^^^, cf. Eq. (p7|), on the two- 
body space restricted by z. The second term comes from 
the one-effective boson exchange (OEBE) and has a form 
similar to (|l5|). 
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Ptc 



ba be 



the same 



(31) 



except for the form factors fx in vertices and the overall 
limitation of the momenta by z [not indicated explic- 
itly in Eq [dw]. Each of these two terms (i.e. projec- 



tion of H 



A(2j 



and voebe) behave for k,k' m like 



the Coulomb potential ( |17| ) with form factors f\ that 
limit changes of the fermion kinetic energies. One can 

<(|)o|Av^|(^o>/Eo 




and z. Numerical calculation confirms that for z > A 
there is no noticeable z-dependence of this matrix ele- 
ment. Fig. 1^ shows how the matrix ele ment depends on 
\ ior z — Qo. As expected in Section V A for small A, 
there is some A-dependence in the result. It emerges be- 
cause at too small As the similarity factors fx start to 
limit the Hamiltonian in the momentum region that is 
important for the bound state formation, and the deriva- 
tion of H\ cannot be carried out precisely using the per- 
turbative renormalization group procedure down to so 
small As. When A and z are large enough, the correction 
Ai^^^-* approaches a fixed finite value, that depends on a. 
This happens because the wave function has a width 
a — and limits the integration over both momenta in 
the matrix element {(f>o\Avx\(j)o). Thus, as seen already in 
Section IV B , the first-order correction is small for small 



coupling constants due to the fast fall-off of the Coulomb 
wave function at large momenta, independently of the 
details of Av that one obtains in the approaches 1 or 
2. The correction is small even for a divergent potential 
such as a (5-function. 



FIG. 7: First-order correction to the ground state binding en- 
ergy as a function of A for z = oo, which shows the magnitude 
of corrections to the well known Schrodinger equation with a 
Coulomb potential expected in the approach 2 in QFT. 



approximate the Schrodinger equation with this QFT 
potential by the equation with a Coulomb potential 
plus a correction, and one can estimate the size of the 
correction using bound-state perturbation theory. For 
this purpose, the difference between potentials vr\ and 
vcouiomb is denoted by Avx. The first-order correction. 



AE 



(1) 



h\Avx\ 



is a function of the parameters A 



Therefore, one needs to look at the second order of 
the bound-state perturbation theory to check the self- 
consistency of the effective particle picture and to com- 
pare it to the approach 1. To see that the effective theory 
does not exhibit the consistency problems the approach 1 
exhibited in Fig. ^, one can closely follow here the deriva- 
tion of Eq. (|2^), but with the OBE potential vqbe re- 
placed by vjix- Again, one can ask whether there is a 
logarithmically divergent dependence on kmax- 

It turns out that for finite values of A there is no such 
divergent dependence. One can safely take the limit 
kmax — oo, since A itself already cuts off sums over in- 
termediate states in the correction. 



4am 



(fki 

7K 



fk2 

E2 



(/)0 (fcl ) AU|X A (fcl , fc2 ) 75- Awij^ A (fc2 , ^3 ) '/'O (^3 ) 



(32) 



r 



Here, Awj^A is defined similarly to Aw||, but with vqbe 
replaced by vrx- Numerical results for this matrix ele- 
ment for different values of A (and for the cutoffs kmax ~ 
A — > 00), are shown in Fig. ^ 

First of all, the results in Fig. || can be considered a 
good approximation to the whole 2nd-order correction 



only for A ^ m (i.e. in the right part of the figure). If 
A is comparable to m, the similarity factors fx limit the 
potential vrx and the high-low and low-high corners of 
the potential matrix (Fig. ^ are practically eliminated. 
The correction coming from the lar ge m omentum region 
selected in the integration in Eq. (p2) is therefore also 
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FIG. 8: Dependence of the large-relative momentum contri- 
bution to the second order bound-state mass correction on A, 
for the cutoff A sent to oo. When one works with effective 
fermions, cutoffs can be sent to oo for any given value of A, 
which is not available in the approach 1 for bare fermions, as 
shown in FigJ3. 



reduced and the other of the 2"* parts of the whole cor- 
rection can contribute in more significant ways than they 
do for large A. Hence, for small A, the results given in 
Fig. ^ are not necessarily a good approximation of the 
whole second-order energy-correction. 

Secondly, in practical work, one needs to lower A as 
far down as possible, possibly below m. Thus, Fig. ^ 
provides only evidence for the self-consistency of the ef- 
fective fermion dynamics in which the convergence in the 
large-relative momentum region is secured by the pres- 
ence of A, and the original QFT cutoffs can be safely sent 
to infinity. 

Thirdly, despite the reservations made above, one can 
expect that the exact 2nd-order correction from the large 
momentum region to the bound-state mass is of the same 
order of magnitude as the part given in Fig. It is 
clear then that the 2nd-order contribution from the large- 
relative momentum range is at least one order of magni- 
tude smaller than the Ist-order result, presented in Fig. |^. 
Therefore, it is the small and moderate momentum region 
that decides about the size of variations of observables 
versus A. In view of these three comments the following 
farther remarks can be made. 

It is visible in Fig. |^ in the case of a = 0.01 that the 
bound state energy is a very slowly varying function of 
A when the latter grows above 0.2m/. One could call 
this value of a a non-relativistic coupling, because all 
the important dynamics happens among virtual effective 
particles with momenta much smaller than the fermion 
mass. Consequently, if the form factor f\ is wider than 
0.2m f, its presence is hardly seen in the final result. Nev- 
ertheless, its presence remains to be essential for keep- 
ing the 2nd-order correction finite and making the large- 
momentum region of QFT a small correction (Fig. 0) to 



the Schrodinger picture when A — > cx), instead of diverg- 
ing and invalidating that picture, as it was happening in 
the approach 1 (Fig. The system of effective parti- 
cles is then shown to be self-consistently non-relativistic 
and the Schrodinger equation with Coulomb potential is 
a good approximation of QFT. Qualitatively, this is the 
situation encountered in QED with a = 1/137. 

However, as a increases, one has to make A larger to 
achieve a fair A-independence of the total binding en- 
ergy, by reducing the size of the A-dependent corrections. 
Fig. 1^ shows that the first-order corrections can be quite 
large and they depend on A. For a — 0.1, the correc- 
tion has a non-neglig ible value {AE'^^^ = OMEq) even 
for A = mf. For a = 0.3 and A = m/, AE'^^^ = O.UEo. 
And if one wanted to make a reduction of the theory to 
the two-fermion state for a = 0.6, one would have to in- 
clude momenta much larger than the fermion mass (for 
A = m/, AE^^^ = OASEq, and it still strongly depends 
on A). In this case, the QFT effective two-body equa- 
tion is no longer close to the non-relativistic Schrodinger 
equation with the Coulomb potential. Also, strong A- 
dependence suggests that considerable errors can come 
from the limitation to only second-order perturbative cal- 
culation of the similarity rotation. Besides, one should 
investigate whether perturbation theory for inclusion of 
higher effective Fock sectors is acceptable when solving 
the Schrodinger equation. 

Fortunately, as one can see from Fig. || at the same 
time the second-order corrections coming from large mo- 
menta remain small, even for a = 0.6. This is an impor- 
tant result, because it suggests that preparing numerical 
procedures for solving effective particle dynamics with a 
few Fock sectors and seeking better starting points for the 
bound-state perturbation theory than the pure Coulomb 
picture, is a good idea to try. Namely, if the region of 
large relative momenta were significant (and not man- 
ageable in perturbation theory, contributing too much in 
the second order), there would be no reason for why lim- 
iting the number of Fock sectors should be a legitimate 
approximation. The uncertainty principle would rather 
suggest that if momenta larger than m/ are important 
then multiparticle states are also important. Our result 
says, that thanks to f\ one can expect the effective Fock 
space expansion to be a legitimate strategy in the ap- 
proach 2, although it requires checking. 

Finally, the difference between the initial cutoff A and 
the parameter A is important for a self-consistent inter- 
pretation of the theory. The cutoff A is an artificial pa- 
rameter that chops off the high-energy part of the initial 
QFT Hamiltonian. One wants to send A to infinity and 
one chooses counterterms in H'^ in such a way that ob- 
servables are independent of A in the limit A — > oo. The 
same huge cutoff appears in the Tamm-Dancoff proce- 
dure with bare particles. There is then no alternative in 
the approach 1 to sending A to infinity in the reduction 
to two-bare fermions Fock sector. This is odd because 
there is no good reason for only a few bare particles to be 
relevant in the bound-state dynamics with the huge mo- 
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mentum cutoff. Indeed, new overlapping divergences are 
obtained this way. Without tools to handle the problem, 
no self-consistent relativistic treatment of bound states of 
fermionic constit uents in QFT can be reached that way 
(see also Section IVA). 

In the conceptually different approach 2, A is absent in 
the effective particle dynamics, and the new finite width 
parameter A is introduced. One adjusts A to contain the 
actual physics of the theory in most economic way, and 
the key measure of success in this respect, besides ver- 
ifying symmetries, is provided by the size of changes in 
the physical results when A is changed. Comparing the 
first and second-order results from Figs. |^ and^, one can 
see that there is a whole range of As and coupling con- 
stants of considerable interest where stability of results 
versus A can be expected in more advanced calculations 
than used in the present pilot study. The renormaliza- 
tion group flow of H\ produces important corrections 
to the well known non-relativistic picture, visible in the 
significant first-order corrections, while the large-relative 
momentum region is self-consistently protected from di- 
vergences. This way, a relativistic treatment of bound 
states of fermions in the Fock space has a chance to be 
developed M. 



VI. CONCLUSIONS 

One can attempt to calculate wave functions of bound 
states of fermions in the Fock space representation of the 
Yukawa QFT in various ways that are described and com- 
pared in the previous Sections on the example of bound 
states of two fermions. Two basic concepts are distin- 
guished. One approach starts from the sector of two 
bare fermions (approach 1), and another one from two 
effective fermions (approach 2). 

Approach 1 leads to overlapping divergences in 
the light-front Hamiltonian dynamics and lacks self- 
consistency in handling the large-relative momentum re- 
gion when one attempts to send the bare cutoff to infinity 
without including infinitely many bare particles. There 
exists an option of removing this defect through sector 
dependent counterterms, but the required construction 
of the full renormalization group triangle with growing 
numbers of bare particles appears overwhelmingly com- 
plex and certainly not completely understood. The basic 
ultraviolet problem comes from short distances in the 
transverse directions and no practical tool exists yet for 
handling huge numbers of bare particles with precision 
required by rotational, parity, and other symmetries of 
the initial Lagrangian. 

Approach 2 is free from the difficulties with large- 
relative momentum convergence of the approach 1. The 
decisive convergence factor is introduced through solv- 
ing renormalization group equations for effective parti- 
cles. The solution includes form factors of width A in 
the interaction vertices and the form factors suppress 
the large momentum domain. This is verified in nu- 



merical estimates described in previous Sections. The 
well-known one-boson-exchange (OBE) potentials that 
are deduced from the on-shell S-matrix elements, are re- 
placed by new one-effective-boson-exchange (OEBE) po- 
tentials, and additional interactions, that are derived in 
the Hamiltonian form independently of the correspond- 
ing S'-matrix. Moreover, one can take advantage of sub- 
sequent S'-matrix calculation in choosing free finite parts 
of the counterterms. Then the renormalized Hamiltoni- 
ans in the QFT Fock space, include relativistic and off- 
energy-shell corrections in accordance with first princi- 
ples of quantum mechanics, special relativity, and renor- 
malization group. However, these principles are satisfied 
only approximately, since the solutions for H\ are found 
only order by order in perturbation theory, and diago- 
nalization of the reduced Hamiltonians Hn is carried out 
numerically. The approximate character is not a weak- 
ness, however, because corrections can be systematically 
investigated and reduced. 

The accuracy of the approximate treatment can be 
estimated by inspecting variation of results when one 
changes the renormalization group parameter A. Exact 
theory would exhibit no dependence on A. In the ap- 
proximate treatment, variations depend on the size of 
the coupling constants. 

For a smaller than about 0.3, there is a wide range of As 
in the Yukawa theory with massless bosons, in which the 
results for the two-fermion bound state mass are stable 
and do not differ much from results obtained from the 
nonrelativistic Schrodinger equation with Coulomb po- 
tential. These values of a can be called non-relativistic. 
For a larger than about 0.3, one has to allow A to 
grow up to the size on the order of my, to achieve A- 
independence of the corrections to the bound-state mass. 
This means that relativistic momenta do matter and the 
non-relativistic Schrodinger equation is not a good ap- 
proximation of the effective dynamics. Still, the effective 
theory is well contained in the range determined by A 
and one can look for solutions of the eigenvalue prob- 
lem without making non-relativistic approximation. For 
the first-order bound-state perturbation-theory correc- 
tions to the non-relativistic approximation are quite con- 
siderable, while the large-relative momentum region in 
second-order corrections contributes merely about 10%. 
Being limited by A, it does not introduce any diverging 
contributions and one does not and should not attempt 
sending A to infinity, in contrast to the cutoff A ^ oo in 
the approach 1. 

All effective-particle calculations are carried out here 
in the light- front form of Hamiltonian dynamics, preserv- 
ing its kinematical covariance under boosts. This form 
enables one to separate the relativistic motion of bound 
states from their internal constituent dynamics. Thanks 
to this separation, one can reduce the description of the 
binding mechanism to the Schrodinger equation for in- 
ternal motion of the constituents without choosing a spe- 
cific frame of reference. All eigenvalue equations derived 
here are independent of the total momentum of the states 
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they describe. Thus, the boost symmetry allows one to 
understand moving bound states in arbitrary motion as 
soon as one understands them at rest, which can help in 
understanding their structure in the infinite momentum 
frame, the opportunity not available in the standard form 
of Hamiltonian dynamics. In fact, the effective particle 
dynamics is in principle well-prepared to tackle Poincare 
symmetry issues in the Fock space including both the 
kinematical and dynamical transformations ^ . 

When the effective light-front theory is defined in a 
perturbative way, the vacuum structure remains simple, 
i.e. b\\0) = d\\0) = a\\0) = 0. In the calculations re- 
ported in this paper, in Yukawa theory, no vacuum struc- 
ture appeared. In the case of theories where spontaneous 
symmetry breaking may occur, one has to include new 
terms in the initial Hamiltonian that contain the infor- 
mation about the possibility that some symmetry may be 
broken spontaneously . No such terms were introduced 
here. 

It is certainly possible to apply the reduction scheme 
described here to other theories than Yukawa. Especially 
interesting case would be QCD, reduced to a nonrelativis- 
tic Schrodinger equation for bound states of constituent 
quarks, since such physical picture is the basis of phe- 
nomenology of hadrons in the particle data tables. It 
is not precisely clear how such picture could follow from 
the canonical QCD itself in any other approach. The fact 
found here that in the Yukawa theory one could enhance 
the coupling constant even up to 0.3 or 0.6 and the con- 
stituent picture remained convergent in the large-relative 
momentum region, is encouraging, because in QCD one 
expects the running coupling to be on the order of 1 for 
widths approaching 1 GeV, cf. ^. 

The discussion of bound states of effective fermions 
would not be complete without answering the question 
of how the binding of effective particles such as nucle- 
ons might emerge from a theory such as QCD. When 
the latter is already written in terms of effective quarks 
and gluons with small A, the Hamiltonian matrix in the 
effective-particle basis can be cut into block pieces that 
lie along the diagonal, and touch each other with their 
corners (see Fig. |9|), some matrix elements being left out. 
Each of the small blocks can be treated as a separate the- 
ory and solved. As a result, one obtains eigenstates or 
scattering states which are closest to true solutions when 
their energies lie somewhere in the middle of a block. The 
accuracy of the solutions diminishes towards the corners 
where the small triangle parts should provide increasingly 
important perturbations from couplings to the neighbor- 
ing blocks. Fig. ||is drawn in the basis of effective quarks 
and gluons. Baryons and mesons in QCD would cor- 
respond to the approximate solutions in the individual 
blocks. It is natural to call the blocks channels and the 
calculation including the interactions between the blocks 
a coupled-channel theory. The transition from the origi- 
nal degrees of freedom (here, effective quarks and gluons) 
to the reduced degrees of freedom (baryons and mesons) 
is achieved by choosing the new basis of states: those 




FIG. 9; Qualitative representation of the idea of solving the 
band diagonal problem in effective Fock state basis in a se- 
quence of diagonalizations of "small" matrices. The width of 
the shaded areas qualitatively represents energy range of the 
couplings between basis states obtained from a new diagonal- 
ization step. 



given by the solution to the channels dynamics when they 
are treated as separate. Thus, the lowest mass channel 
with baryon number 2 contains only nucleons interact- 
ing through potentials, the next channel with invariant 
masses between 2mN + TOtt and 2mN + 2mT^ may include 
two nucleons and one pion, etc. The bound or scatter- 
ing states from higher blocks could be approximated by 
products of the eigenstates from lower-mass blocks. From 
the mathematical point of view, the transition step is 
a change of basis in describing eigenstates of relatively 
small matrices (in comparison to the initial Hamiltonian 
matrix) and taking into account the residual interactions 
between different channels later to refine the basis states. 
The refinement step would be executed in the new basis 
using the same type of division into blocks, but with the 
new blocks centered where the previous blocks were least 
accurate. This sequence of steps can be repeated many 
times. Thus, the whole transition process to the effective 
theory written in terms of new degrees of freedom would 
form a sequence of patterns indicated in Fig. |^. The 
shaded shapes in Fig. ^ indicate the boundaries of non- 
zero matrix elements of the effective Hamiltonian ma- 
trices between the incoming and outgoing states, which 
qualitatively indicates their range in energy. 

By the way, it is also natural to see in similar terms 
the transition from electrons and nuclei to atoms, or 
from atoms to molecules. The universal character of this 
scheme is the change of basis in the Hilbert space of states 
in quantum mechanics, and the Hamiltonian scheme for 
effective particles does open this option to detailed stud- 
ies in all these theories. Once the width is sufficiently 
small, it slows down, or freezes degrees of freedom that 
require larger mass changes than A to be excited. One 
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can also imagine that reducing the width farther in the 
new basis, a next level of effective degrees of freedom can 
be isolated using the same scheme. If the reduction were 
started with a theory underlying QCD, i.e. at extremely 
large A, the scope of the reduction scheme would be even 
greater. 

As a final remark let us indicate that quantum field 
theories such as Yukawa theory, or sigma-models, may be 
helpful in discovering the acceptable structure of channel 
Hamiltonians even without making the transition from 
QCD to nucleons and mesons explicitly. This option fol- 
lows from the fact that when the width in the effective 
theory is made small enough, the meson-fermion vertex 
contains a narrow form factor and the emission of mesons 
is weak and slow, even for quite large coupling constants. 
The structure of the small-width effective theory may 
be then rather simple and universal in the sense that 
many initial complex theories may reduce to the same 
small-width effective one. Thus, the apparently academic 
Yukawa theory, or other QFTs of that kind, may become 
very useful in indicating what kind of effective reduced 
Hamiltonians one can expect from the more elaborate 
analysis that derives nucleon-pion dynamics from QCD. 
In fact, a small number of free parameters, such as the 
form factor width, the size of the coupling constant, or 
the size of effective masses in a basic set of terms, may 
be sufficient to reproduce the bulk of the effective model 
predictions. It is also possible that no significant width 
dependence is observed in the model results if the param- 
eters are correlated as if the effective model were derived 
from a Yukawa theory or a sigma-model using only a per- 
turbative expansion for artificially small couplings that 
only later are extrapolated to the required large values 



one can look for R and Hn in perturbation theory in 
Hj. This leads to the lowest (second order) expression 
for Hr: 



Hr = PHP + ^PHiQHiP - ^PlhQHiP 



(A4) 



where the underlining denotes free energy denominators, 
displayed in Eq.(^. 



APPENDIX B: NOTATION FOR SPINORS 

The four-component spinors of fermions and anti- 
fermions in light-front Hamiltonians have the following 
form. 



1 



^mpX 



1 



[p+A+ + A_(m + a-Lp^)] vx ,(B1) 



where the projection operators are given by 
the spinors of particles of zero velocity equal 



(B2) 




(B3) 
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APPENDIX A: REDUCTION PROCEDURE 

The eigenvalue equation for a Hamiltonian H, 

Hm = (Al) 

for some low eigenvalues, is replaced using the operation 
R pi], E2l by an eigenvalue equation for eigenstates \(p) — 

P + R^ R\tp) of the reduced Hamiltonian Hr, given by 
the following formula, 

' ' - ' - ' (A2) 



Hr = 



{p + r^)h{p + r:}- 



R'^R 



Vi + R^iR VT 

If one splits the initial Hamiltonian into the free and 
interaction parts, 



H^Hn + Hr 



(A3) 



and the two-component spinors are 
XI - ( J ) , X-i 



(B4) 



The matrices appearing on the right hand side of Eq. 
(Bl) represent Lorentz transformations belonging to the 
small group that preserves the four-vector rj defining the 
light front = through the condition rjx — , where 
X is an arbitrary four- vector. 

Fermion field operator fulfilling free Euler-Lagrange 
equations ipm{x) is expanded for light-front "time" = 
in terms of creation and annihilation operators as: 



per ^mpcr^ 



(B5) 



APPENDIX C: LIGHT-FRONT MOMENTA 

For two particles of momenta pi and p2 , the total mo- 
mentum P and the relative momenta x, k"*" are defined 
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as follows. 



or, equivalently, 



Pt 


= xP+ , 


(CI) 


1 


Pt 


= il-x)P+ , 


(C2) 




pi 


= XP^ + , 


(C3) 




pk 


= (1 - X)P^ ~ . 


(C4) 





1 + 



(CIO) 



Unless stated otherwise, 



(C5) 



The shorthand notation in the integrals involves the mea- 
sure 



2(27r)3p+ 



(C6) 



Creation and annihilation operators are normalized by 
the commutation relation 



ap,4 



2(27r)3/c+j3(fc_p) 



(C7) 



for bosons, and in analogous way by anticommutation 
relations for each kind of fermions. 

In order to exhibit relationship between QFT and the 
well-known equal-time Schrodinger equation, one uses 
relative momentum k 129(1, 



k^ = 



x(l — x) 



(C8) 
(C9) 



APPENDIX D: SCHRODINGER'S SOLUTION IN 
COULOMB POTENTIAL 



Equation 



(the Schrodinger equation for positronium without spin 
and with reduced mass = m/2) has the ground state 
eigenvalue 



R 1 2 



and the normalized ground state wave function 

1 



(a2 + F)2 ' 
with N = ^Jia^ jjF' /it and a — afi. 



(D2) 



(D3) 
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